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Abstract 

Every route of a transport network approaching equilibrium can be represented by a vector of Euclidean 
space which length quantifies its segregation from the rest of the graph. We have empirically observed that 
the distribution of lengths over the edge connectivity in many transport networks exhibits scaling invariance 
phenomenon. We give an example of the canal network of Veneice to demonstrate our result. The method 
is applicable to any transport network. 
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Transport networks are used to model the flow of 
commodity, information, viruses, opinions, or traf- 
fic. They typically represent the networks of roads, 
streets, pipes, aqueducts, power lines, or nearly any 
structure which permits either vehicular movement or 
fiow of some commodity, products, goods or service. 
The major aim of the analysis is to determine the 
structure and properties of transport networks that 
are important for the emergence of complex flow pat- 
terns of vehicles (or pedestrians) through the network 
such as the Braess paradox [T] . This counter- intuitive 
phenomenon occurs when adding more resources to a 
transportation network (say, a new road or a bridge) 
worsens the quality of traffic by creating longer delays 
for the drivers, rather than alleviate it. The Braess 
paradox has been observed in the street traffic of New 
York City and Stuttgart, [1]. 

In the present Letter, we show that while approach- 
ing equilibrium, a transport network can be embed- 
ded into Euclidean space M^"-'^, N being a number 
of vertices. Then, every edge of the network is rep- 
resented by a vector which length quantifies its seg- 



regation from the rest of the graph. We have empir- 
ically observed that the distribution of lengths over 
the edge connectivity in urban transport networks ex- 
hibits scaling invariance phenomenon. The relation 
between the connectivity of city spaces and their cen- 
trality known as intelligibility is a key determinant of 
human behaviors in urban environments, [3]. 

In most of researches devoted to the improving of 
transport networks, a primary graph representation 
of urban networks is used in which streets and routes 
are considered as edges of a planar graph, while the 
traffic end points and street junctions are treated as 
nodes. The usual city map based on Euclidean ge- 
ometry can be considered as an example of primary 
city graphs. 

However, another graph representation can be use- 
ful if we are interested in describing the transport net- 
work at equilibrium. Given a connected undirected 
graph G{V,E), in which V is the set of nodes and E 
is the set of edges, we introduce the traffic function 
/ : E ^ (0,00 [ through every edge e G E. It then 
follows from the Perron-Frobenius theorem 4 that 



the linear equation 

/(e) - E fie')exp{-hi{e')), (1) 

where the sum is taken over all edges e' € E which 
have a common node with e, has a unique positive 
solution /(e) > 0, for every edge e ^ E, for a fixed 
positive constant h > and a chosen set of positive 
metric length distances £{e) > 0. This solution is nat- 
urally identified with the traffic equilibrium state of 
the transport network defined on G, in which the 
permeability of edges depends upon their lengths. 
The parameter h is called the volume entropy of the 
graph G, while the volume of G is defined as the sum 
Vol(G) — \ X)ee_E^(e)- The volume entropy h is de- 
fined to be the exponential growth of the balls in a 
universal covering tree of G with the lifted metric, 

The degree of a node i ^ V is, the number of its 
neighbors in G, degQ(z) = ki. It has been shown in [8] 
that among all undirected connected graphs of nor- 
malized volume, Vol(G) = 1, which are not cycles and 
for which ki ^ 1 for all nodes, the minimal value of 
the volume entropy, min(/i) = i X^iGV ^* ^"-"S ~ 1) 
is attained for the length distances 
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where ki and kj are the degrees of the nodes linked 
hy e G E. It is then obvious that substituting ([2]) and 
min(/i) into ([1]) the operator exp {—M{e')) is given by 
a symmetric Markov transition operator, 
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where i and j are the nodes linked by e' S E, and 
the sum in ^ is taken over all edges e' € E which 
share a node with e. The symmetric operator ([3]) 
rather describes time reversible random walks over 
edges than over nodes. In other words, we are invited 
to consider random walks described by the symmetric 
operator defined on the dual graph G*. 

The Markov process ([3]) represents the conserva- 
tion of the traffic volume through the transport net- 
work, while other solutions of ([T]) (with h > mm{h)) 



are related to the possible termination of travels 
along edges. If we denote the number of neighbor 
edges the edge e £ E has in the dual graph G* as 
= deg(3*(e), then the simple substitution shows 
that w(e) = y/q^ defines an eigenvector of the sym- 
metric Markov transition operator defined over the 
edges E with eigenvalue 1. This eigenvector is pos- 
itive and being properly normalized determines the 
relative traffic volume through e G i5 at equilibrium. 
Eq.Q relates the equilibrium transport flows on the 
graph G to the stationary distribution of random 
walks defined on its dual counterpart G* and empha- 
sizes that the degrees of nodes are a key determinant 
of the transport networks properties. 

The notion of traffic equilibrium had been intro- 
duced by J.G. Wardrop in [3] and then generalized 
in [10] to a fundamental concept of network equilib- 
rium. Wardrop's traffic equilibrium is strongly tied 
to the human apprehension of space since it is re- 
quired that all travellers have enough knowledge of 
the transport network they use. Dual city graphs are 
extensively investigated within the concept of space 
syntax, a theory developed in the late 1970s, that 
seeks to reveal the effect of spatial configurations on 
the human perception of places and behavior in urban 
environments, [HI [3]. Spatial perception that shapes 
peoples understanding of how a place is organized de- 
termines eventually the pattern of local movement, 
which is quantified by the space syntax measure be- 
ing nothing else, but an element of a transition prob- 
ability matrix of a Markov chain [12j , with surprising 
accuracy jlBI . Random walks embed connected undi- 
rected graphs into the Euclidean space K''^"^. This 
embedding can be used in order to compare nodes 
with respect to the quality of paths they provide for 
random walkers and to construct the optimal coarse- 
graining representations. 

While analyzing a graph, whether it is primary 
or dual, we assign the absolute scores to all nodes 
based on their properties with respect to a trans- 
port process defined on that. Indeed, the nodes of 
G{V, E) can be weighted with respect to any mea- 
sure m = ^i^ymili, specified by a set of posi- 
tive numbers > 0. The space i'^{m) of square- 
assumable functions with respect to the measure m 
is the Hilbert space Ti-iV). Among all measures which 
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can be defined on V, the set of normalized measures 
(or densities), 1 = X^igy""*!*; of essential inter- 
est since they express the conservation of a quantity, 
and therefore may be relevant to a physical process. 

The fundamental physical process defined on a 
graph is generated by the subset of its linear auto- 
morphisms which share the property of probability 
conservation; it can be naturally interpreted as ran- 
dom walks. The linear automorphisms of the graph 
are specified by the symmetric group S^v including 
all admissible permutations p G Sat taking i V to 
p{i) € V and preserving all of its structure. 

Markov's operators on Hilbert space form the nat- 
ural language of transport networks theory. Being 
defined on a connected aperiodic graph, the transi- 
tion matrix of random walks is a real positive 
stochastic matrix, and therefore, in accordance to the 
Perron- Frobenius theorem [4] , its maximal eigenvalue 
is 1, and it is simple. The left eigenvector ttT =^ tt 
associated with the eigenvalue 1 is interpreted as a 
unique equilibrium state tt (the stationary distribu- 
tion of random walks). The Markov operator T is 
self-adjoint with respect to the normalized measure 
associated to the stationary distribution of random 
walks TT, 



f =1(^1/2^^-1/2 



TT ' 1 TT ' 



(4) 



where is the transposed operator. In the theory 
of random walks defined on graphs [14] and in spec- 
tral graph theory [TS] , basic properties of graphs are 
studied in connection with the eigenvalues and eigen- 
vectors of sclf-adjoint operators defined on them. 
The orthonormal ordered set of real eigenvectors ipi , 
i = 1 ... TV, of the symmetric operator T defines a 
basis in H(F). 

In particular, the symmetric transition operator T 
of the random walk defined on connected undirected 
graphs is Tij — 1/ y^kikj, if z ~ j, and zero other- 
wise. Its first eigenvector tpi belonging to the largest 
eigenvalue /ii = 1, 



V'lT = -01, V'm 



(5) 



describes the local property of nodes (connectivity), 
since the stationary distribution of random walks is 



T^i = ki/2M where 2M — J^iev^i- remaining 
eigenvectors, {ips}'^^2^ belonging to the eigenvalues 
1 > M2 ^ • ■ ■ IJ-N ^ — 1 describe the global connect- 
edness of the graph. For example, the eigenvector 
corresponding to the second eigenvalue fj,2 is used to 
define the spectral bisection of graphs; it is called the 
Fiedler vector if related to the Laplacian matrix of a 
graph [T5] . 

Markov's symmetric transition operator T defines 
a projection of any density a € 'H{V) on the eigen- 
vector -01 of the stationary distribution tt. 



crT = xjji + a± T, 



(6) 



in which a± is the vector belonging to the orthogonal 
complement of ■01. In space syntax, we are interested 
in a comparison between the densities with respect 
to random walks defined on the graph G. Since all 
components "01. ^ > 0, it is convenient to rescale the 
density a by dividing its components by the compo- 
nents of ipi , 

^^ = ^- (7) 

Thus, it is clear that any two rescaled densities 
CT,p G 'H{V) differ with respect to random walks 
only by their dynamical components, (a — p)T* = 
—p±)T*', for all t > 0. Therefore, we can de- 
fine the distance || . . . ||t between any two densities 
established by random walks by 
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or, using the spectral representation of T, 

II ||2 >^ (ct_L - P± |0s)(0s I CT± - PJ_) 
h-p\\T = 2^ Y— ' 

(9) 

where we have used Diracs bra-ket notations espe- 
cially convenient for working with inner products and 
rank-one operators in Hilbert space. 

If we introduce a new inner product for densities 
CT,p^n{V) by 
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then © is nothing else but || cr — p | 
||p||^-2((T,p)j,, where 
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is the square ol the norm ol G T^iV) with respect 
to random walks defined on the graph G. 

We finish the description of the Euclidean space 
structure of G induced by random walks by mention- 
ing that given two densities cr, p € Ti-iV), the angle 
between them can be introduced in the standard way, 



cos Z {p, cr) 



\P\\t 



(12) 



The cosine of an angle calculated in accordance to 
(fT2| has the structure of Pearson's coefficient of lin- 
ear correlations. The notion of angle between any 
two nodes of the graph arises naturally as soon as we 
become interested in the strength and direction of a 
linear relationship between the flows of random walks 
moving through them. If the cosine of an angle (fT2)) 
is 1 (zero angles), there is an increasing linear rela- 
tionship between the flows of random walks through 
both nodes. Otherwise, if it is close to -1 (tt angle), 
there is a decreasing linear relationship. The correla- 
tion is {tt/2 angle) if the variables are linearly inde- 
pendent. It is important to mention that as usual the 
correlation between nodes does not necessary imply a 
direct causal relationship (an immediate connection) 
between them. 

In order to illustrate the approach, we have stud- 
ied five different patterns of compact urban trans- 
port networks. Two of them are situated on islands: 
the street network in Manhattan and the network of 
Venetian canals. We have also considered two me- 
dieval German cities developed within the fortresses: 
Rothenburg ob der Tauber in Bavaria and the down- 
town of Bielefeld in Eastern Westphalia. To supple- 
ment the study of urban canal networks, we have also 
examined it in the city of Amsterdam. In all analyzed 
transport networks, we observe that the segregation 
of a node measured by ifTTj) scales negatively with its 
connectivity; the slopes of the regression lines slightly 
exceed 2. 



In the present Letter, we describe the Euclidean 
space structure of space syntax for the spatial net- 
work of 96 Venetian canals which serve the function 
of roads in the ancient city that stretches across 122 
small islands. While identifying a canal over the plu- 
rality of water routes on the city map of Venice, the 
canal-named approach has been used, in which two 
different arcs of the city canal network were assigned 
to the same identification number provided they have 
the same name. 

In accordance to (fTTt . the density 1^, which equals 
1 at i €z V and zero otherwise, acquires the norm 



associated to random walks. Its square, 
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expresses the access time to a target node in random 
walk theory [14j quantifying the expected number of 
steps required for a random walker to reach the node 
i ^ V starting from an arbitrary node chosen ran- 
domly among all other nodes with respect to the sta- 
tionary distribution tt. 

The notion of spatial segregation acquires a statis- 
tical interpretation with respect to random walks by 
means of ([T3l) . In urban spatial networks encoded by 
their dual graphs, the access times (fT3|) strongly vary 
from one open space to another and could be very 
large for statistically segregated spaces. It is remark- 
able that the norm a canal of Venice acquires with 
respect to random walks scales with its connectivity 
(see Fig. HI). 

The Euclidean distance between 1^ and Ij induced 
by random walks. 
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is the commute time in theory of random walks being 
equal to the expected number of steps required for a 
random walker starting at z G ^ to visit j € V and 
then to return back to i, [T^ . 

Indeed, the structure of vector space induced 
by random walks cannot be represented visually, how- 
ever if we choose a node of the graph as a point of 
reference, we can draw the 2-dimensional projection 
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Figure 1: The scatter plot of the connectivity vs. the 
norm a node in the dual graph representation of 96 Vene- 
tian canals acquires with respect to random walks. Three 
data points characterized by the shortest access times rep- 
resent the main water routes of Venice: the Lagoon of 
Venice, the Giudecca canal, and the Grand canal. Four 
data points of the worst accessibility are for the canal sub- 
network of Venetian Ghetto. The slope of the regression 
line equals 2.07. 

of Euclidean space by arranging other nodes at the 
distances and under the angles they are with respect 
to the chosen reference node. The 2-dimensional pro- 
jection of the Euclidean space of Venetian canals set 
up by random walks drawn for the the Grand Canal of 
Venice (the point (0,0)) is shown in Fig. [2] Nodes of 
the dual graph representation of the canal network in 
Venice are shown by disks with radiuses taken equal 
to the degrees of the nodes. All distances between 
the chosen origin and other nodes of the graph have 
been calculated in accordance to ([14]) and p2)l has 
been used in order to compute angles between nodes. 
Canals positively correlated with the Grand Canal of 
Venice are set under negative angles (below the hor- 
izontal), and under positive angles (above the hori- 
zontal) if otherwise. 

The radiuses of disks display the equilibrium con- 
figuration of flows along the Venetian canals when 



Figure 2: The 2-dimensional projection of space syntax 
of Venetian canals built from the perspective of the Grand 
canal of Venice chosen as the origin. The labels of the 
horizontal axes display the expected number of random 
walk steps. The labels of the vertical axes show the degree 
of nodes (radiuses of the disks). 



the trafRc volume is conserved. It is evident from 
Fig. [5] that disks of smaller radiuses demonstrate a 
clear tendency to be located far away from the origin 
being characterized by the excessively long commute 
times with the reference point (the Grand canal of 
Venice) , while the large disks which stand in Fig. [2] 
for the main water routes are settled in the closest 
proximity to the origin that intends an immediate 
access to them. 

Probably, the most important conclusion of space 
syntax theory is that the adequate level of the pos- 
itive relationship between the connectivity of city 
spaces and their integration property (vs. segre- 
gation) called intelligibility encourages peoples way- 
finding abilities [3l . Intelligibility of Venetian canal 
network reveals itself quantitatively in the scaling of 
the norms of nodes with connectivity shown in Fig. [1] 
and qualitatively in the tendency of smaller disks to 
be located on the outskirts of the Venetian space syn- 
tax displayed in FigO 
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